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Abstract. We study the direct and an inverse source problem for the radiative transfer equation 
arising in optical molecular imaging. We show that for generic absorption and scattering coefficients, 
the direct problem is well-posed and the inverse one is uniquely solvable, with a stability estimate. 



1. Introduction 

We consider an inverse source problem arising in optical molecular imaging (OMI) which is 
currently undergoing a rapid expansion. The design of new biochemical markers that can detect 
faulty genes and other molecular processes allows us to detect diseases before macroscopic symptoms 
appear. This has been studied extensively in the bioengineering literature. See for instance [3], 
[5], [7]. Unlike higher-energetic markers used in classical nuclear imaging techniques such as single 
photon emission computed tomography (SPECT), positron emission tomography (PET) as well as 
magnetic resonance imaging (MRI), optical markers emit relatively low- frequency photons. The 
objective of OMI is to reconstruct the concentration of such markers from their radiations measured 
at the boundary of the domain. The radiations in OMI are governed by the equations of radiative 
transfer and the inverse problem in OMI is thus an inverse transport source problem, at least once 
the optical properties of the underlying medium are known. We now describe more precisely the 
mathematical problem. 

We assume that O is a bounded domain of R n with smooth boundary. We will assume also that 
is strictly convex. This is not an essential assumption since for the problem that we study, one 
can always push the boundary away and make it strictly convex, without losing generality. 
The radiative transport equation 

(1) 9- V x u(x,0) + a(x,9)u(x,9) - I k(x,9, 9')u(x,9') dO' = f(x), u\ d _ sn = 0. 

where the absorption a and the collision kernel k are functions with a regularity that will be 
specified below. The source term / is assumed to depend on x only. 

In section [2] we study the direct problem. We show that for an open and dense set of absorption 
and scattering coefficients the direct problem ([1]) is well-posed. See Theorem [1] for details. 

The boundary measurements are modeled by 

(2) Xf(x,9) = u\ d+S n, (x,0)ed+Sn 

where u(x,9) is a solution of (pQ), and (9+5T2 denotes the points x £ dfl with direction 9 pointing 
outwards. 

In section 3 we consider the inverse source problem, that consists in determining the source term 
/ from measuring Xf. Notice that in the case a = k = the linear operator X is the standard 
X-ray transform and when k = 0, X is a weighted X-ray transform (see section [2]). 
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This inverse problem has been considered in several papers in the mathematical and engineering 
community [2], [8], [TT], [13], [TT] . In particular in [2] it is shown that one can prove uniqueness 
when k = k(x, 0-9'), and k is small enough in a suitable norm. We show that for the absorption and 
scattering in an dense and open subset we can uniquely determine the source / from the boundary 
measurements. We also prove a stability estimate. See Theorem [2] for details. 

2. The Direct Problem 

Set 

(3) T = 9-V x , Ti=T + a, T = T + a-K, 

where a is viewed as the operator of multiplication by o~(x,9), and K is the integral operator in 
Let u solve 

(4) Tu = f, u\ d „ sn = 0. 

As mentioned in the introduction the operator X is the X-ray transform if a and k are both zero 

,0 

Xf(x,9) = If(x,6) := / f(x + t6)dt, (x,6)ed+Sn (<r = fc = 0), 

Jr-(xfi) 

where ±r±(x, 9) > are defined by (x, x + t±(x, 9)) € d±<ST&. We will always extend / as outside 
f2, therefore we can assume that we integrate above over R. If k = 0, then X reduces to the 
following weighted X-ray transform 

(5) Xf(x,9) = I a f(x,9) := J E(x + t9,9)f(x + t9)dt, (x,9)ed + Sn (k = 0), 
where 



(6) E(x, 9) = exp f- y o-(x + se,f)da 

If o" > depends on x only, this is known as the attenuated X-ray transform, that is injective, and 
there is an explicit inversion formula, see [TO], pQ. 

We define the adjoint X* of X w.r.t. the measure dE defined above. We will view X as a 
perturbation of and our goal is to show that X*X is a relatively compact perturbation of I^I a . 

First we will analyze the direct problem. In the next theorem, / is allowed to depend on 9 as 
well. 

Theorem 1. There exists an open and dense set of pairs (cr,k) £ C 2 (£l x S' 11 ^ 1 ) x C 2 (£l x 5 n_1 x 
S 1 ™ -1 ), including a neighborhood of (0,0), so that for each (o~,k) in that set, 

(a) the direct problem ^ has a unique solution u E L 2 (f2 x S 11 ^ 1 ) for any f € L 2 (£l x S" 1-1 ) 
depending both on x and 9. 

(b) X extends to a bounded operator 

X : L 2 (n x S"" 1 ) — ► L 2 (d + Sn, dS). 

Proof. We start with the analysis of the direct problem In what follows, let To, T\ and T 
denote the operators given by (HP in L 2 (Q x S 1 ™" 1 ) with domain 

D(T ) = £)(Ti) = D(T) = {/ G L 2 (n x S^ 1 ); 9-V x u£ L 2 (n x S^ 1 ), n| a _ Sf , = 0} . 
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We will assume here that / depends both on x and 9. Note first that the solution to the problem 
([4]) with k = is given by u = T^ 1 f, where 

(7) [T- 1 f](x,6) = J exp^- J a(x + T9,9)dT^j f(x + s9,9)ds. 

This follows easily from the fact that E is an integrating factor, i.e., To = E~ l T\E. 
Apply Tf 1 to both sides of (gj) to get 

u = T 1 - 1 (Ku + f). 
We therefore see that (|3J) is equivalent to the integral equation 

(8) (Id - T^K)u = Tf 1 /- 

Therefore, if Id — T^ 1 K is invertible, ([!} is uniquely solvable, and the solution is given by 

(9) u = T- 1 f = (ld-T- 1 Kr 1 T- 1 f. 
We also get that T _1 exists if and only if (Id — T^ l K)~ l exists. Then 

(10) u = T~ l f = (Id - T{ l K)- l T^ l f. 
When / depends on x only, set 

(11) [Jf}(x,9):=f(x). 
Then 

(12) u = T~ x Jf = (Id - Tf 1 if)" 1 Tf 1 Jf. 
Lemma 1. The operator KT± J : L 2 (fi) — > L 2 ($7 x S 1 ™ -1 ) is compact. 
Proof. Let first / depend both on x and 9. Then 



(13) 



where 



[KT^f}(x,9) = [ k{x,9,9') I expf- f° a{x + t9' ,9') dr\ f(x + s9',9')dsd9' 

J S™- 1 J-oo \ Js J 



\x — y\ n 1 \~ \x — y 



f [ y, t- — -r ) d y> 



E(x,s,e') = exp (- J a(x + T9',9')dT^j 



(we replaced s by — s and then made the change x — s9' = y). 

Note that when / depends on 9 as well there are no enough integrations above to show that 
T{~ K is compact. On the other hand, if / depends on x only, i.e., if we have Jf above with such 
an /, we have 

(14) [KT± J]f(x, 9) = / " yl n _^ l ^f(y)dy. 

Jn \ x — y\ 

The integral above is a typical singular operator with a weakly singular kernel, and an additional 
parameter 9, see [9], [2]. Under the smoothness assumptions on a, k, it is easy to see that 
OqKT^ 1 and dxKT^ 1 are bounded operators, see Proposition [1] below. This completes the proof 
of the lemma. □ 
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As we mentioned above, the arguments above do not prove that KT X is compact on L (Q x 
5 n_1 ). On the other hand, its square is compact. 

Lemma 2. The operator T^KT^ 1 : L 2 (Q x S 71 ' 1 ) -» L 2 (Q x S 11 ' 1 ) is compact. 
Proof. Replace /(y, rfzff) in ([TSJ by 

[if/] (v, ) = / fc f y, f^r , ^ /(v, d(9'. 
V \ x -y\J Js"' 1 V F-y| / 

Then the compactness follows from the same arguments as in Lemma [TJ Indeed, we have 

x-y 



r r a l x, y, \x — y\, i _ . , 
[^r i - 1 ET/](x,0)= // . " ' f(y,8')dydd' 

iisixs"- 1 \ x - y\ 

with an obvious definition of a. In particular, all second order derivatives of a are bounded. 
Let g(x,9,9') be the y-integral above, i.e., the r.h.s. above becomes J g(x, 9, 9') d9'. Then by 
Proposition [TJ below, 

r \d x g(x,9,9')\ 2 dx<C [ \f(x,9')\ 2 dx 



for any 9, 9'. In particular, 

\d x g(x,9,9')\ 2 dxd9' <C\\f\\ 2 L2 . 



II 



inxS"- 1 
Then 

2 



dxd9 



d x KT{ 1 Kff= I! [ d x g(x,9,9')d9' 

<C ff [ \d x g(x,9,9')\ 2 d9'dxd9 



<C'\\f\\ 2 



2 

La- 
It is easy to see that dgKT^ 1 Kf E L 2 as well. This, and the estimate above, imply the compactness 
of KT{ l K. □ 

We proceed with the proof of part (a) of the theorem. We are looking for k so that T _1 exists. 
Consider 

A(X) = (id - (XKT{ 

in L 2 (£l x S 71 ^ 1 ). The operator (i^T^ -1 ) 2 is compact, and for A = 0, the resolvent above exists. By 
the analytic Fredholm theorem [12], A(X) is a meromorphic family of bounded operators. In partic- 
ular, it exists for all but a discrete set of A's. Thus for the those A's, the resolvent (id — XKT^ 1 ) 
exists and is given by 

(15) (Id - XKT^y 1 = (Id + XKT{ 1 ) A (A). 

Indeed, it is obvious that the operator on the r.h.s. above is a right inverse to Id - XKT^ 1 . For 
|A| <C 1, one can use Neumann series to show that it is left inverse as well. One can use analytic 
continuation around the poles to show that this remains true for all A that are not poles. 

By ([9]), then T -1 exists for such A's and k replaced by Xk. In particular, this shows that the 
set of such (k, a) is dense. Standard perturbation arguments show that the set of k : s for which 
Id — XKT^ 1 is invertible, is open in C° for a fixed a; and the set of pairs (a, k) € C° x C° with the 
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same property is open, too. Since we just showed that it is dense as well in C° x C°, this completes 
the proof of (a). 

We proceed with the proof of (b). For X we get, see ([9]), 

(16) Xf = R+T- l f = R+(Id - T^Ky^f, 
where 

R+h = h\ d+sn . 

If / depends on x only, then 

(17) Xf = R+T-\jf = R+(Id - T{ l K)- l T^\jf. 
Notice first that 

(18) (Id - T^Ky 1 ^ 1 = T-f x (Id - KT^y 1 , 

and in particular, the resolvent on the left exists if and only if the resolvent in the r.h.s. does. We 
therefore have 

Xf = R + Ty(ld-KTyy 1 f. 
To prove (b), it is enough to show that 

R+ry : L 2 yi x s n - v ) — ► L 2 (d + Sfl, d£) 
is bounded. A straightforward computation (see also [I]) shows that 

f [ f(x-td,e)dt6S= [ f(x,6)dx 

for any / G L 1 (fi x S"™ -1 ). Therefore, 



de 



\R+Tyf\\ 2 LHd sn ^= [ \R + Tyf(x,9)\ 2 dX< [ f f(x + t8,8)dt 

Ja+sn Jd+sn Jr-(x,e) 

<[ (\t-(x,6)\ f \f(x + t9,9)\ 2 dt ] dS 
Jd+sn \ Jt-(x,o) I 



2 

dS 



'8+Sfl V Jt-(x,i 
< diam(^) ||/|||2 (nx5 n-i)- 
This completes the proof of Theorem [TJ □ 

3. The Inverse Source Problem 

In this section we consider the inverse source problem. The next theorem shows that for generic 
(a, k) there is uniqueness and stability. As mentioned in the introduction a similar result has been 
proven in [2] in the case where k = k(x, 0-9'), and k is small enough in a suitable norm. 

Fix a bounded domain f^i so that Q\ 3> fi. Extend (a, k) with regularity as below to functions 
in S7i with the same regularity. We chose and fix that extension as a continuous operator in those 
spaces. Define the operator X\ : L 2 (yL\) — > L 2 (d+SMi) in the same way as X. We will be interested 
in the restriction of X\ to functions / supported in £1. We always extend such / as zero to Q\. 
This corresponds to taking measurements on dQ\ instead of dQ. 
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Theorem 2. There exists an open and dense set of pairs 

(19) (a, k) € C 2 (ti x S" 1 - 1 ) x C 2 (ti x x S™- 1 ; C 71 * 1 ^" 1 )) , 

including a neighborhood of (0,0), so that for each (o~,k) in that set, the conclusions of Theorem^ 
hold in Oi, and 

(a) the map X\ is injective on L 2 (£l), 

(b) the following stability estimate holds 

(20) \\f\\v m <C\\X*Xf\\ mm , V/GL 2 (tt), 
with a constant C > locally uniform in (cr,k). 

Remark 1. The smoothness requirement on k can be reduced to k E C 2 if k is of a special form, 
like k = Q(9)k(x, 6') or a finite sum of such, see ([26]) . ([27|) . 

^From now on, we will drop the subscript 1, and all operators below are as defined before but in 
the domain Vt\. We assume that (a, k) are such that T -1 exists. We assume now that X is applied 
to / that depends on x only. For now, it is not important that / is supported in Q; that will be 
needed in (|3Tj) and after that; so we apply X to functions in L 2 (Qi). By (fT71) . 

(21) X = I a + L, L := R + (-Id + (Id -T^Ky 1 )^ 1 J, 
see also ([5]). Then 

(22) X*X = r a I a + £, £:= I*L + L*I a + L*L. 

In our analysis, we will apply a parametrix of I%l a to X*X. That parametrix is a first order 
operator. For this reason, we study d x I*L. 

Lemma 3. The operators 

d x I*L, 8 x L*I a , d x L*L 
are compact as operators mapping L 2 (£li) into L 2 (Qi). 

Proof. To analyze I^L, note that L also admits the following representation 

(23) L = R+T^KT^ 1 (id - KTf 1 ) ~* J. 

We need to study I*R + T^ 1 KT^ 1 h, where h = h(x, 9). Notice first that 

(24) [I*h](x)= [ E(x,9)h$(x,9)d9, 



where E denotes complex conjugate, and h* is the extension of h € C{dj r SQ.\) as a constant along 
the lines originating from x in the direction —9, see e.g., [6J sec. 4]. In other words, h^(x,9) = 
h(x + t + (x,9),9). Next, R+T± h looks just like I a , see ([5]) but with / there depending on 9 as 
well. Therefore, 



\I*R + T{ l \g{x) = [ E(x,9) [ E(x + t9,9)g(x + t9,9) dt 
This yields (see [6] again): 

[i:R + T 1 - 1 g](x)= 11 E(x,9)(Eg)(x + t9,9)d9dt 



d9. 



(25) 



|y — x\ n 1 
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where F even (x, 9) is the even part of F as a function of 9. If we set g = KT X h, that would give us 
I*R+T{ l KT{ l h. 

Instead of assuming (|19p . we will make the following weaker assumption at this point: k can be 
written as the infinite sum 

oo 

(26) k(x,9,9') = ^Q j {9)K j {x,9') 
with some functions Qj and Kj so that 

oo 

(27) ^lie.n 

3=1 

One such way to do this is to choose Qj to be the spherical harmonics Yj\ then Kj are the correspond- 
ing Fourier coefficients. Then \\Yj H^i^n-i) < C(l + \j), where X 2 are the eigenvalues of the positive 
Laplacian on S" 1 " 1 . Since Xj = 0(j 1 ^ n ~ 1 ^), for the uniform convergence of (|26[) it is enough to have 
||«ilU- < C(l + Aj)- n - £ , e > 0. This would be guaranteed if k E L°° (fij x S™- 1 ; C^+V™ -1 )) 
by standard integration by parts arguments. Therefore, the hypothesis (|19p of the theorem implies 

HMD, j27D. 

Under this assumption, for KjT^~ 1 h, where i£j has kennel QjKj, we have, see f 1 1 3 [) . 
[K J T 1 - 1 h](x,9) = e j (9)[B j h](x), 

F-yr V \x-y\J 

We claim now that Bj (Id - KTf x ) J : L 2 (fti) -> L 2 (^i) is compact. We have 

(Id - KT- 1 )' 1 J = J + (Id - KT^y 1 KT~ X J. 

By Lemma [H the second term on the right is compact. Therefore, it remains to show that BjJ is 
compact. Observe that BjJh is given by (|28p with h = h{x). The compactness then follows from 
Proposition [TJ assuming that Kj € C 2 . On the other hand, BjJ is compact under the assumption 
that Kj € L°° only, by O Theorem VII. 3. 3]. Moreover, its norm is bounded by C||/Cj||l«>. 
We can write now 

d x P a L = d x r a R + T{ l KTy (Id - xryy 1 J 

(29) 00 r _ii 

= J2 [d^R+T^QjJ] [Bj (Id - KTf 1 ) J . 

3=1 

We notice first that d x I*R + Ty@jJ : L 2 (fti) -► L 2 (^i) is bounded by Proposition QJb), compare 
to ([25]) . with a norm bounded by C 1 1 c" 1 1 c 2 1 1 ®* j 1 1 J^" 1 - The operator in the second brackets is compact 
by the claim above. Therefore, each summand in the r.h.s. of (|29p is a compact operator with 
a norm not exceeding C||0j ||#i \\Kj \\l°°, where C depends on a as well. Then the series in (|29p 
converges uniformly by (|27p . Under this condition, d x I*L is compact. 



s 
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To analyze d x L*L, we will follow the proof above. It is enough to show that d x L*R + T 1 1 ®jJ '■ 
L 2 (0,) — » L 2 (Q,i) is bounded. We have, see 



d x L*R + T~ 1 Q j J = d x (R + T~ l (id - KT~ l ) XT' 1 J) R+T^QjJ 

( 3 °) ' */ — 1 \ * 

= 8 X (KT^J) (ii+Tf 1 (Id - KTf 1 ) ) R+T^&jJ. 



Since R+T-^ 1 is bounded, it remains to show that the operator d x (KT^ x J)* : L 2 (f2i x S n l ) — > 
L 2 (Q) is bounded, as well. The kernel of the latter is, see (|14p . 

(x,( y ,e)) —> d x ^ 1 ly { , v 

\y — x\ n 1 

Then the boundedness of d x (-fTT^ 1 ,/)* then follows as in Lemma [2j 

Finally, the fact that d x L*I a is bounded follows from the proof for d x L*L. Indeed, d x L*I a = 
d x L* R + T^ 1 J ', compare with ([30]) . where we can set Qj = 1. 

This completes the proof of Lemma [3j □ 

Proof of Theorem^ We return to the analysis of the operator X*X, see (|22p . We showed in 
Lemma[3]that, up to a relative compact operator, X*X coincides with I*I a . Assume that a and k 
are C°°. Let Q be a parametrix (of order 1) to the elliptic ^DO I* la in Qi . We restrict the image 
of Q to L 2 (f2), i.e., we view Q as an operator Q : — > L 2 (f2). Then for any / supported in 
f2, we have 

(31) QI*U = f + K ± f, 
where K\ is of order —1 near f2. Apply Q to X*X to get 

(32) QX*Xf = f + K 2 /, K 2 := + Q£. 

Then K2 : L 2 (J7) — > L 2 (0) is compact. We get that the problem of inverting X* X is reduced to a 
Fredholm equation. We will show, that it is generically solvable, as in the theorem. 
We show first that the set of pairs for which X is injective is dense. 

By the results in [6], see Theorems 1 and 2 there, if a is real analytic in a then I a is 
injective, and therefore is injective as well. Moreover, a small C 2 (£l), perturbation preserves 
that property. Actually, the remark after the [BJ Theorem 2] shows that this is true even for 
small enough C 1 perturbations. Fix a real analytic in Fix k as well so that (a, k) belongs 
to the generic set in Theorem dj related to fij, and the regularity assumption (fT9|) is satisfied. 
That can be done for an open dense set of /c's by the proof of Theorem [TJ Consider X related 
to (<r, Afe) with A belonging to some complex neighborhood C of [0,1]. The operator K2 in (|32p 
depends meromorphically on A € C. Indeed, K\ is related to (a, k) (i.e., to A = 1), and is therefore 
independent of A. The parametrix Q is also independent of A. The analysis above shows that C is 
a meromorphic function of A because L has that property, see (|15p and (|2ip . For A = 0, we have 
C = 0, and then K2 = K\. By adding a finite rank operator to Q, we can arrange that Id + K\, 
see (j3T|) . is injective, see also the proof of Proposition 4 in [15]. We can then apply the analytic 
Fredholm theorem again in C with the poles of (Id — XK)^ 1 ^ 1 removed. The latter is a connected 
set, containing A = and A = 1. The analytic Fredholm theorem then implies that QX*X is 
invertible for all A in that set with the possible exception of a discrete set. In particular, there are 
A's as close to A = 1 as needed with that property. For those A's, X*X and X are injective as well. 
This shows that there is a dense set of pairs (a, k) in the space f|19|) so that X is injective. Lets us 
call that set U. 
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We show next that for (a, k) in some neighborhood of U, X is still injective. 

Let (k,a) G U. Then X : L 2 (n) -» L 2 (<9^i,d£) is injective . Then X*X : L 2 (0) -> JET 1 ^) is 
injective as well, as an integration by parts shows. By adding a finite rank operator to Q, we can 
arrange that Id + K±, see (|3ip . is injective, as above. Then Id + K\ is invertible on L 2 (Q), and we 
deduce that ([H holds. 

The analysis above implies that the norm ||^*-X"||i / 2(n)-+fl"i(n 1 ) depends continuously on (a, k) as 
in (fl9l) . Therefore, we can perturb (a,k), and (|20|) would remain true because the perturbation of 
the r.h.s. will be absorbed by the l.h.s. On the other hand, injectivity of X*X implies injectivity 
of X. 

This proves that the set of pairs (o~,k), for which X is injective, is open subset of the (generic 
set) of pairs, for which the direct problem is guaranteed to be uniquely solvable by Theorem [TJ 
Moreover, (120p with C locally uniform. 

This completes the proof of Theorem [21 □ 

In the proof of the theorem, we used the following proposition about singular operators. 
Proposition 1. Let A be the operator 



/a [x, y, \x — y\ 
V | g _ y |»-i ' ' y /(y)dy 

with a(x,y,r,8) compactly supported in x,y. Then 

(a) If a £ C 2 , then A : L 2 — > H 1 is continuous with a norm not exceeding C\\(x\\q2. 

(b) Let a(x,y,r,9) = a'(x,y,r,6)4>(6). Then \\A\\ L 2_^ H i < C\\a'\\ C 2\\4>\\ H i( S „-iy 

Proof. We recall some facts about the Calderon-Zygmund theory of singular operators, see [9]. First, 
if K is an integral operator with singular kernel k(x,y) = <f)(x,8)r~ n , where 6 = {x — y)/\x — y\, 
r = \x — y\, and if the "characteristic" (j) has a mean value as a function of 9, for any x, 
then K is a well defined operator on test functions, where the integral has to be understood in the 
principle value sense. Moreover, K extends to a bounded operator to L? with a norm not exceeding 
C sup x \\<t>(x, -)|ll, 2 (S n -i)) see IS Theorem XI. 3.1]. The characteristic (j) does not need to have zero 
mean value in 6 but then the integral has to be considered as a convolution in distribution sense. 
The latter is well defined because the Fourier transform of the kernel w.r.t. the variable z = r6 is 
homogeneous of order 0, thus bounded. 

Also, see Theorem XI.11.1], if B is an operator with a weakly singular kernel ip(x, 0)r~ n+1 , 
then d x B is an integral operator with singular kernel d x [(3(x, 6)r~ n+1 ]. The latter, up to a weakly 
singular operator, has a singular kernel of the type 4>r~ n , and the integration is again understood 
in the principle value sense, see the next paragraph. In particular, the zero mean value condition 
is automatically satisfied. 

In our case, (3 = a depends on y and r as well. Assume first that it does not, i.e., B is as above. 
Extend (3 as a homogeneous function of 9 of order near S 11 ^ 1 . Then 

gM) _ n J iRj ^d(3/de j de j /3 Xi (x,e) 



Xl r n— 1 v ' r n^ / j r n—l Q x . 

j 1 

-^- n )^l j + 2^ r n ^ ~ W + r n-l 



r 

3 



(i-n)e l f3 + df3/de i p Xi (x,o) 
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We used the fact that ^ ■ 9jdf3/d9j = because (3 is homogeneous of order in 9. It is not hard to 

show that the "characteristic" <p(x, 9) = {l—n)9i(3+d(3 / d9i has zero mean over S^ -1 , see [9j p. 243]. 
In this particular case (a(x,y,9) = f3(x,9), independent of y, r), statement (a) can be proven as 
follows. Choose a finite atlas of charts for S"™ -1 so that for each chart, n—1 of the 9 coordinates (that 
we keep fixed in R") can be chosen as local coordinates. By rearranging the x, and respectively, 
the 9 coordinates, in each fixed chart, we can assume that they are 9' = (9i, . . . , 9 n —i)- Then 
d(3/d9 n = — Y17=l 9f3/d9i. Then in (f33l) . we have derivatives of j3 w.r.t. 9' (and x) with smooth 
coefficients. The contribution of the first term then can be estimated by the Calderon-Zygmund 
theorem. The second term is a kernel of a weakly singular operator. The following criterion can be 
applied to it: If K has an integral kernel k(x, y) with the property 

(34) sup / \k(x,y)\dx < M, sup / \k(x, y)\dy < M, 

x J y J 

then K is bounded in I? with a norm not exceeding M |16|, Prop. A. 5.1]. 
This proves (a) for a = j3. 
To replace (3(x,9) above by a(x,y,9), write 

a(x, y, r, 9) = a(x, x, 0, 9) + rj(x, y, r, 9). 

To prove (b), write first as above, 

(35) a(x, y, r, 9) = p'(x, 9)<j){0) + ry(x, y, r, 9)<j)(9), (3'(x, 9) := ai(x, x, 0, 9), 
where 7 G C . Notice then that in (|33p . with (3 = (3'cb, we have 

(1 - n)9if3 + 313/ d6i = (1 - n)9iP'(f) + Qdp'/dOi + f3'd(f>/d9 i . 

Choosing local coordinates as above, and applying the Calderon-Zygmund theorem again, we get 
that the first term above contributes a singular operator with a norm not exceeding ||ai||c*i H^Hh 1 - 
The second term generates an operator with a kernel j(x,y,r,9)(j)(9)r~ n+ ' 2 . Differentiate w.r.t. 
x, and we still get a weakly singular operator which norm can be estimated as in (|34p to give a 
norm not exceeding 1 1 71 1 1 c i \ \ cj) \ \ H i . □ 

Remark 2. The only second order derivatives of a that were needed in the proof of (a) were 
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